Proton structure corrections to hyperfine splitting in muonic hydrogen 
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We present the derivation of the formulas for the proton structure-dependent terms in the hyper- 
fine sphtting of muonic hydrogen. We use compatible conventions throughout the calculations to 
derive a consistent set of formulas that reconcile differences between our results and some specific 
terms in earlier work. Convention conversion corrections are explicitly presented, which reduce the 
calculated hyperfine splitting by about 46 ppm. We also note that using only modern fits to the 
proton elastic form factors gives a smaller than historical spread of Zemach radii and leads to a 
reduced uncertainty in the hyperfine splitting. Additionally, hyperfine splittings have an impact on 
the muonic hydrogen Lamb shift /proton radius measurement, however the correction we advocate 
has a small effect there. 



I. INTRODUCTION 

The recent conundrum with the proton charge ra- 
dius [l| inspires reconsideration of the corrections that 
enter into determinations of proton size dependent energy 
sphttings in hydrogen atoms, particularly muonic hydro- 
gen. This note is about the determination of the polariz- 
ability and elastic corrections to the hyperfine splitting in 
muonic hydrogen, where there is a disagreement implicit 
in the literature which should be resolved. The numer- 
ical consequences of the resolution have a small effect 
upon current measurements, but could be consequential 
when accurate measurements of muonic hydrogen hyper- 
fine splitting become available and if other sources of 
uncertainty are reduced. 

For the leading order proton structure-dependent part 
of the hyperfine splitting (HFS) in muonic hydrogen, 
the calculation as commonly implemented requires, we 
claim, further correction due to overlapping terms be- 
tween the elastic and polarizability parts of the calcula- 
tion. The situation occurs because terms can be moved 
from the inelastic to the elastic part of the calculation us- 
ing the Gerasimov-Drell-Hearn-Hosada-Yamamoto sum 
rule (which relates the proton anomalous moment to in- 
elastic cross sections) [2rS|i ^-^id details of the move arc 
handled differently in different sources. For cases known 
to us, the differences are of no numerical consequence 
for ordinary hydrogen, but are noticeable for muonic hy- 
drogen. Articles on this subject often calculate just the 
elastic or the just inelastic corrections 0, and some- 
times both in the same article [1, |^ . A calculation of the 
proton structure corrections that combines elastic and 
inelastic (under the heading of "polarizability") results 
from incompatible sources will obtain a deficient total, 
and such calculations appear to exist in practice. At a 
practical level, given that many full muonic HFS articles 
{e.g., [13, Eli) s-lso usefully compile higher-order correc- 
tions and small terms from other effects, perhaps the best 
procedure is not to start from scratch, but to add to the 
existing calculations the needed term that coordinates 
the elastic plus inelastic total. 



Part of the current motivation for discussing the HFS is 
the conundrum occasioned by the recent measurement of 
the proton charge radius using the Lamb shift in muonic 
hydrogen The HFS has some, albeit small, effect 
there. The Lamb shift is the energy difference, mostly 
due to QED effects, between the 25*1/2 and 2Pi/2 hydro- 
gen levels. More precisely, given that each of the lev- 
els is hyperfine split, it is the energy difference between 
the properly weighted average of the hyperfine levels of 
the 25*1/2 and 2Pi/2 states. The goal is to use the mea- 
sured Lamb shift to infer the extra overall shift of the 
25*1/2 energy due to proton size corrections. The exper- 
iment actually measures the energy difference between 
the 25*^^^^ state (where F is the total lepton -I- proton 
angular momentum) and the IP^p^ states, and relies 
upon calculation of the 25i/2 HFS and the 2P3/2-2P1/2 
splitting to obtain the Lamb shift. The P-states involve 
no proton structure corrections to the accuracy required. 
The 25f/=i state is shifted up by (1/4) of the full HFS, 
so a change in the HFS has an effect upon the Lamb shift 
extraction. However for the correction discussed in this 
paper, the effect translates into less than 1% of the cur- 
rent discrepancy between the muonic hydrogen value of 
the charge radius [l| and the CODATA value which 
is based mainly on the ordinary hydrogen Lamb shift and 
confirmed by electron scattering data [T3| . 

Our outline for this paper is as follows. Many of the 
formulas we use were given in our earlier work Q, 
but without derivations being shown. Not much detail 
wasgiven as well for the muonic polarizability calculation 
in [7[. Since there are some differences to be assessed, we 
show the calculation in moderate detail in Sec. [Hi It 
parallels the well-known calculation for negligible lepton 

mass EMI, 

but now keeps the mass general. A reader 
who does not wish to check the details on a first reading 
can proceed directly to Sec. IHIl where we isolate the 
term that should be added to some existing muonic HFS 
calculations, evaluate it numerically, and compare its size 
to other known terms. The effect of this change to the 
HFS upon the current measurement of the proton charge 
radius from the Lamb shift is shown in Sec. IIVI A short 
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discussion is offered in Sec. |V] 



and the spin-dependent part of Wui, is 



II. HYPERFINE SPLITTING 

The relevant diagram is shown in Fig.m 
We neglect the Fermi momentum of the atomic lepton, 
and note that the proton leg plus photons is the same 
as off-shell forward Compton scattering on the proton, 
which is calculated from the matrix element of a time 
ordered product 
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A e^''-HpS\TjAaiAO)\pS), (1) 



where is the electromagnetic current and the states 
are normalized by {pS\p'S) = 2E{2'k)^ 5^ {p — p') (for the 
same 5") . We need only the spin dependent part of T^i, , 
which is antisymmetric in its indices and is expanded 
using the scalar functions Hi and H2, 



(2) 



p-q 



- S-q p^ 



Qa = ^Z*^ is given in the proton rest frame, qo — p ■ q/rup, 
and is the proton spin vector. 
From the definitions 



Hi{-qo,Q^) 
H2i-qo,Q^) 



Hi{qo,Q'), 
-H2{qo,Q^). 



(3) 



The functions Hi 2{qo,Q^) are not known either from 
ab initio calculation or measurement. However, the func- 
tions can be constructed from their imaginary parts. The 
imaginary parts come from the contribution where the 
intermediate electron and hadron states, in Fig. [1] are 
on-shell physical states, which in turn are squares of 
electron-proton clastic or inelastic scattering amplitudes. 
Hence we can obtain -ffi,2('Z05 Q^) using dispersion rela- 
tions and data on electron-proton scattering cross sec- 
tions. 

In terms of standard notation, electron-proton scatter- 
ing cross sections are given from the hadronic tensor 

= ^ /^'^ {ps\ MO,M0)] \ps) , (4) 



p. 




FIG. 1: Diagram giving proton structure-dependent correc- 
tions to tlie hyperfine splitting. 



Hence (for qo > 0), 

ImHi(go,Q') 
Imi/2(go,g') 



(5) 



q p'^ 



52(90,(9^) 



— .9l('?0,Q^) , 
10 



(6) 



The spin-dependent structure functions gi and 172 are 
measured (over some kinematic range) in inelastic po- 
larized e-p scattering at laboratories including CERN, 
SLAC, HERMES, ELSA, and JLab. For elastic scatter- 
ing they are given in terms of the elastic form factors 
by 
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elastic 
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elastic 



= -lrF2GMSil-x) , 



(7) 



where x — Q^/{2p ■ q), Fi{Q^) and F2{Q^) are the Dirac 
and Pauli form factors, respectively, and Gm and (for 
completeness) Ge are 



Gm — Fi + F2 , 

Ge ^ Fi — -^-^^2 
4to^ 



(8) 



The functions Hi^ilOiQ^) have poles in \qQ\ at the 
elastic point and cuts beginning at the inelastic threshold 



i^th = m-TT + {ml + (9^)/(277ip) 



(9) 



and going to infinity along the real axis. We use the 
integration contour illustrated in the complex q^ plane 
in Fig. [5J assume zero contribution from the part of the 
contour at infinite distance, and obtain the dispersion 
relations 



Hiiqo,Q^) 



2mp q^F,iQ^)GMiQ^) 
TT (g2 + ie)2 _ Am^q^ 



dv 



2 9i{v,Q ) , 



(10) 



H2{qo,Q^ 



2mpmpiyF2iQ^)GMiQ^) 
TT (g2 + ig)2 _ Am'^qQ 

2mpqo f°° dv 1 / ^2n /..n 



Evaluating the two-photon, structure-dependent, con- 
tributions to the HFS, Fig. [l] one obtains after making 
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where the elhpses stand for functions that may depend 
on and other variables but not on the individual com- 
ponents of Q. The auxiliary functions [IE [Hi are 



-3r + + 2(2 - r)Vr(T + l), 



/32(t) = l + 2r-2V^(7TT): 
and have limits 



/32(t) 



f+0(l/r), r^oo, 



1 + T^O, 
+ C'(l/r), T^oo. 



(16) 



(17) 



FIG. 2: Contour in complex gg plane for Cauchy integral. 



a Wick rotation, 
-£^2°^ arag ( d^Q 1 

~E7 " (1 + Kp)^2 y 02 Q4 + 4^2g2X 

X ((2g2 + gg)i/i(*go,g^)- sg^g^^^^^^^'^'^ 



(12) 



where is the Fermi energy 
i'V^ = — MBMp - 



3 MB (1 + m^/mp)^ ' 



(13) 



and TTir is the lepton reduced mass. The Wick rotation 
entails the notations 



When the lepton mass can be neglected inside the 
integrals, the results (jTS)) are directly applicable with 
T = IQ^ . When the lepton mass cannot be neglected, 
one has denominators that split into partial fractions. 
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(g4 + 4m2g2)(g4 + 4m2g2) 

1 / 1 
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(g4 + 4m2g2)(^.2+g2) 

1 



1 



1 



2(g2-4m2T) V^Q' + Qo nQ"^ 



4g4(m2 - Tn\) \TpQ^ + Ql t,Q^ + Ql) ' 
1 



(18) 
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q 



«go , 

i dQo d^q = i d'^Q . 
2 _ ^2 



q 



(14) 



Hence the results including lepton mass are given in 
terms of the same auxiliary functions and the notations 



After substituting for -ffi,2 using the dispersion rela- 
tions, one has a five dimensional integral. The structure 
functions 51,2(^^7 g^) depend only on two variables, and 
three integrals can be done analytically. The results are 



d^Q 2g2 + Ql 
d^Q Ql 



dg2 1 

^/52(r)...,(15) 



4m2 



4mi 



(19) 



are useful. 

The contribution of the two-photon diagram to the 
HFS becomes 
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rpbox 



amirrip 



Ep 2(1 + Kp)7r(m2 - nif) 



(/32(Tp)-/32(T,)) F2{Q^)GMiQ') 



2amf 



(1 + Kp)7r 



dQ^ 



dv 



)2 _ 4jrijT 



Mr) 



dQ^ 



dv 



12 /^^^ Q2 _ Am^.T 



5(/32(r)-/32(r,)) 52(i^,Q') 



(20) 



Conventionally, the proton structure-dependent eor- 
rections are split into three terms, 



As = + Afl + Apoi = 



Tpbox o 



^ f21) 



called Zemach, recoil, and polarizability terms. 

We have subtracted from the box diagram the iteration 
of the lowest order one-photon exchange diagram, since 
in a bound state calculation that contribution is already 
included [Tl| . This cancels the infrared divergence in 
the box diagram. The visible effect of the subtraction is 
to give the "—1" in the Zemach term displayed below. 
Including this subtraction gives the final result, up to a 
reorganization of terms. 

The Zemach term is 



8am,. 



dQ (Ge{Q'')Gm{Q^) 



-2amrrz , 



(22) 



where rz is the Zemach radius. The first part of the 
Zemach term is obtained from the first line of Eq. ([20)) 
with /3i (Ti ) replaced by the first term in its low argument 
limit and Fi replaced by G^;. The Zemach term is finite 
in the nonrelativistic limit, meaning the proton mass go- 
ing to infinity with electron mass and proton size held 
fixed. 

Next continuing to the gi term, notice that aside from 
the overall mi factor, the integral diverges in the 
limit. To demonstrate this, note that the dv integral 
for — in this limit is just a numerical factor times 
the left hand side of the Gerasimov-Drell-Hearn-Hosada- 
Yamamoto sum rule Q-Q, 



Anir 



dv 2 
— gi(v,0) = -Kp. 



(23) 



The dQ^ integral then diverges at its lower end. 

Conventionally, this near divergence is removed by 
adding an extra term to the gi integral. The gi and 
g2 terms together go into the polarizability correction, 
and with standard notation. 



ipol 



ami 



2(1 + Kp)nmp 



(Ai + A2) 



(24) 



with 



PiirdF^iQ') 



dv Q^P^{T)~Amjv^p^{n) 2^ 

+ 4m„ / — — — ^^-^ gi{v,Q ) 



"j.^.v^ Q^-Amjv 



Ay = -12™: 



°° dQ^ 

""Jo 

V^ g4 _ 4„-j2^2 



g2{v,Q^). (25) 



With the F2 term, Ai is finite in any limit. The rea- 
son for multiplying F2 by /3i(r^) is to make it compati- 
ble with [11,0]. Other choices could be made, including 
just using the numerical factor (9/4) as in the purely 
electron case [H, [l^ , but every choice has consequences 
elsewhere. 

What remains is the recoil term, 



T^rn^ Jo 



I + Kp 



amimp 



2(1 + Kp)7r(m2 - mj) ^ Jo 



•p 

°° ^r)2 



n 



FiiQ^)GM{Q^) 



dQ' 



ami 



(/32(rp)-/32(T,)) F2{Q')Gm{Q^) 



2(1 + Kp)-Kmp Jq 



2 /3i(t,)F|(Q2). (26) 



The first term comes because the Zemach correction was 
written in terms of Ge instead of Fi . The middle terms 
can be recognized in Eq. (|20p . The last term compensates 
the extra term in the polarizability correction. 

Let us now recall a well-known albeit old-fashioned way 
to calculate the elastic contributions alone. That is to 
calculate the box and crossed box diagrams. Fig. [31 using 
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a photon-proton vertex given by 

r.(g) = 7.^1 (Q') + 7^(y.pq'F2{Q 



(27) 



for incoming photon momentum q. We do not advocate 
doing the calculation this way because the above vertex 
cannot be complete or correct when the intermediate pro- 
ton is off-shell. One can correctly obtain the imaginary 
part of the elastic part of the Compton amplitude this 
way, since only on-shell configurations give the imaginary 
part. However, this calculation has been done for the full 
box, and the results are well-known and widely used 
{e.g. [13, Eli) for the elastic terms. 

The dispersive calculation, barring questions of sub- 
tractions, is complete and correct, and using the term 
described above ensures that the elastic terms — defined 
as the Zemach plus recoil terms — are the same as the 
result of calculating the box diagrams directly with the 
vertex of Eq. (P7)) . To state the same thing in other 
words, if one uses the elastic contributions to the HFS 
found in [s], Q one must use the polarizability terms given 
here in order to obtain the correct total proton structure 
contributions to the HFS. 



III. CONVENTION CONVERSION 
CORRECTIONS 



As noted, the mi ^ result for the polarizability cor- 
rection in [7| differs from ours. The difference is entirely 
in the F| terms. They have 



with 



l3o{Te)F2{Q'^) + (other terms same) 



/3o(t) =2vM7TT)-2t. 



(28) 



(29) 



This result is correct if used to complement the elastic 
terms calculated in the conventions of, for example, [l7| . 

We can, however, find two examples in the literature 
where the Ref. Q polarizability corrections are combined 
with the elastic terms of Refs. [1, Q (the Az + as 
displayed here) [HI ■ In this case one should add a 
further correction for muonic hydrogen HFS given by 



Apoi(corr.) 



2(1 -t- Kp)7rm 





(30) 



^71 




P p p p 

FIG. 3: Elastic box and crossed box. 



which evaluates to —46 ppm, or about —1.0 fieV. To 
two figures (in the ppm number), the result is the same 
using the form factors from Arrington, Melnichouk, and 
Tjon [il], Kelly [l^^rrington and Sick or the new 
Mainz form factors mHH. 



IV. EFFECT UPON THE LAMB SHIFT AND 
PROTON RADIUS 



In the context of the Lamb shift, the proton structure 
correction to the HFS that we discuss in this paper is 
small compared to the larger proton structure correction 
of the 25 energy level and has only a slight effect on the 
extracted proton radius. 

The muon Lamb shift experiment measures the transi- 
tion between the 25^^'^^ and 2P^^^ states, and following 
the discussion in the introduction one may write an equa- 
tion from Ref. [ij as 

E{2P^,=^) - £'(25'f/=i) = {rp independent terms) 



- 5.2262 meV- 



fm 



0.0347 meV- 



fm^ 



2S 
HFS 



(31) 



where rp is the proton charge radius and the "r^ inde- 
pendent terms" are well calculated. 

A simple 46 ppm reduction in the 2S'-state HFS results 
in an increase in rp, but one that is smaller than the un- 
certainty limit of the experiment [l[ . Stated differently, if 
one wants to make the above equality work with the CO- 
DATA value for the radius, one needs to find a 310 /leV 
or so extra contribution to the energy, which we do not 
have from this source. 



V. CLOSING REMARKS 

We have reconsidered the hyperfine splitting for 
muonic hydrogen, advocating a unified calculation in or- 
der to clearly track how each term is defined. We re- 
ported numerical results earlier, and have here detailed 
the derivation and shown an additional term that should 
be included if the 0{a^) elastic and polarizability terms 
are taken from specific different publications. 

The result of this term is to reduce the quoted hyper- 
fine splitting by about 1.0 ^eV. 

As a remark regarding the totality of the proton struc- 
ture dependent corrections, many modern form factor fits 
give a Zemach radius toward the larger end of its former 
range [H, _24 2^ , and with a smaller spread [l^] . Since 
the spread of values for the Zemach radius contributed 
the largest uncertainty to the calculated HFS, there could 
be a notable reduction in the quoted error limit. The 
first three form factors listed in Table U show this pos- 
sibility. Table U was prepared using our own results for 
the 0{a^) (lowest non-trivial order) proton structure cor- 
rections and using the extensive compilation of [llj for 
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the QED, higher-order, and additional small corrections, 
leaving these other terms in [Tl| untouched and quoting 
results for the 2S'-state for definiteness. We have also 
indicated the Zemach radius that follows from each form 
factor parameterization. 

TABLE I: Hyperfine splitting for the 25 state of muonic hy- 
drogen, using different modern analytic fits in tlie terms tliat 
involve elastic form factors. 



Form factor fit Ejfps (meV) rz (fm) 

AMT [18] 22.8123 1.080 

Kelly [19| 22.8141 1.069 

AS [20| 22.8105 1.091 

Mainz 2010 [21-23| 22.8187 1.045 



However, the latest electron scattering Mainz re- 
sults [llj create an exception. They agree well with the 
CODATA value [l2| for the charge radius but give a mag- 
netic radius that is noticeably smaller than previous, but 
still modern, fits. Along with this, they give a smaller 
Zemach radius [2^ [2^ than the other form factor fits in 
Table H and increase the spread of calculated values for 

77^25 
^HFS- 

Nearly all the uncertainty in the HFS calculation 
comes from the proton structure terms. Quoting from Q , 
the uncertainty in the polarizability contribution is ± 114 
ppm or ± 2.6 /xeV. Determining the uncertainty in the 
elastic contributions from the spread of results from the 



selection of form factor parameterizations gives ±4.1 ^eV 
(±180 ppm). This would have been ±1.8 /xeV (±80 
ppm) before the latest result. Adding in uncertainties in 
quadrature gives ±4.9 /^eV, meaning that one should be 
allowed already to quote a smaller uncertainty as in 

Ejfps = 22.8146(49) mcV , (32) 

compared with Effp^ = 22.8148(78) meV obtained 
in and used in One may expect that newer 

comprehensive analyses of the magnetic form factor data 
could reduce the uncertainty limit farther. The central 
value quoted is the midpoint of the values in Table U and 
is little moved. 

This note has focused on one correction where we be- 
lieved there were some definite statements to be made, 
and in the light of the PSI experiment on the Lamb 
shift, all discrepancies need to be sorted out. We 
have not reassessed the whole set of corrections to the 
HFS [nl, Hill as has been done recently for the spin- 
independent case [2^ [13] ■ 
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